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Recent progress in the nonperturbativc solution of (3+l)-dimensional Yukawa the- 
ory and quantum electrodynamics (QED) and (l+l)-dimensional super Yang— Mills 
(SYM) theory will be summarized. The work on Yukawa theory has been extended 
to include two-boson contributions to the dressed fermion state and has inspired 
similar work on QED, where Feynman gauge has been found surprisingly conve- 
nient. In both cases, the theories are regulated in the ultraviolet by the inclusion 
of Pauli— Villars particles. For SYM theory, new high-resolution calculations of 
spectra have been used to obtain thermodynamic functions and improved results 
for a stress-energy correlator. 



1. Introduction 

Numerical techniques can be successfully applied to the nonperturbative so- 
lution of field theories quantized on the light cone. 1,2 ' 3 Unlike lattice gauge 
theory, 4 wave functions are computed directly in a Hamiltonian formula- 
tion. The properties of an eigenstate can then be computed relatively easily. 
There have been a number of successes in two-dimensional theories, 3 but 
in three or four dimensions the added difficulty of regulating and renormal- 
izing the theory has until recently limited the success of the approach. 

Here we discuss recent progress with two different yet related approaches 
to regularization. One is the use of Pauli- Villars (PV) regularization 5, 6,7,8,9 
and the other supersymmetry. 10 The particular applications to be discussed 
are to Yukawa theory and QED in 3+1 dimensions with PV fields and to 
super Yang-Mills (SYM) theory in 1+1 dimensions. In the latter case, 
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extension to 2+1 dimensions has already been done; 11 however, the most 
recent developments have used two dimensions as a proving ground. There 
we consider in particular a stress-energy correlator 12 ' 13 ' 14 and analysis of 
finite-temperature effects. 15 

The light-cone coordinates 1 that we use are defined by x^ — x° ± 
x 3 , x± = (x 1 ,^ 2 ), with the expression for x ± divided by \[2 in the case 
of supersymmetric theories. Light-cone three-vectors are denoted by an 
underline: p= (p + ,p±). 

The key elements of the PV approach are the introduction of negative 
metric PV fields to the Lagrangian, with couplings only to null combina- 
tions of PV and physical fields; the use of transverse polar coordinates in the 
Hamiltonian eigenvalue problem; and the introduction of special discretiza- 
tion of this eigenvalue problem rather than the traditional momentum grid 
with equal spacings used in discrete light-cone quantization. 2,3 The choice 
of null combinations for the interactions eliminates instantaneous fermion 
terms from the Hamiltonian and, in the case of QED, permits the use of 
Feynman gauge without inversion of a covariant derivative. The transverse 
polar coordinates allow use of eigenstates of J z and explicit factorization 
from the wave function of the dependence on the polar angle; this reduces 
the effective space dimension and the size of the numerical calculation. The 
special discretization allows the capture of rapidly varying integrands in the 
product of the Hamiltonian and the wave function, which occur for large 
PV masses. 

For supersymmetric theories, the technique used is supersymmetric dis- 
crete light-cone quantization (SDLCQ), 16 ' 10 which is applicable to theories 
with enough supersymmetry to be finite. This method uses the traditional 
DLCQ grid in a way that maintains the supersymmetry exactly within the 
numerical approximation. The symmetry is retained by discretizing the 
supercharge Q~ and computing the discrete Hamiltonian P~ from the su- 
peralgebra anticommutator {Q~ , Q~} = 2\/2P~ . To limit the size of the 
numerical calculation, we work in the large-7V c approximation; however, 
this is not a fundamental limitation of the method. 
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2. Yukawa theory 

The Yukawa action with a PV scalar and a PV fcrmion is 



d 4 x 



(i) 



~\ tyiVfy - (S/^ifr") V>i + miVi^i 



-g(<fo + <h.)Wo + ^1)^0 + ^1)] • 
From this we obtain the light-cone Hamiltonian 

P = E/ ^^^^(-^^tW&MW 



(2) 



+ 



2 — 2 

E / dq^—^(-l) j a](q) aj (q) 

3 

E / d£^{[^-2»(P>£) + lMp + £,«)] b[ s (p)al(q)b l ^ s (p + q) 
[Uj(p,q) + U t (p + q,q)] b[ s {p)a\(q)bi, s {p + q) + h.C.} , 



i,j,k,s 



where antifcrmion terms have been dropped. No instantaneous fcrmion 
terms appear because they are individually independent of the fermion 
mass and cancel between instantaneous physical and PV fermions. The 
vertex functions are 



Uj(p,q) 



fJ 



V2s(p, q) 



9 e 2s • P-L 



(3) 



V167T 3 p + x/q+ v ' ■" \/8iT 3 p + \/q + 

with e*2 S = — -^(2s,i). The nonzero (anti)commutators are 

[a t (q),a}(q^] = (-1)%% - £ '), (4) 

&t s ,( £ ')} = (-1)^,^(5 -pO- 

We construct a dressed fermion state, neglecting pair contributions; it 
takes the form 

<M£) = E 2 A f + (£)|o) + £ / ^.(flXte^teiJIo) (5) 

+E / d ii d %hUg v g 2 ) /T l^ b L(fc)"](fci)4(fc 2 )|o) + • • ■ 
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The wave functions f s {x n , q±n) satisfy the coupled system of equations that 
results from the Hamiltonian eigenvalue problem P + P~<& + = M 2 $ + . Each 
wave function has a total L z eigenvalue of (1) for s = +1/2 (—1/2). 
The coupled equations are 

r p+ 

m*z l + Y / (-±Y +3 P + J dq{f Vj -{q)\y + (P-q,q) + V±(P,q)] (6) 
+ fi>j+ (g) [U V (E -£,£) + Ui (P, £)] } - M 2 Zl , 



m 2 i + q\ + I*} + q\ 



fi 3 s{q) 



i -y y 

+ J2(- 1 f Ui'^AVUE- q,q) + V-(P,q)} 

v 

+Zi> S Si+ [Ui (P -£,£) + U V (P, £)] } 



(7) 



x {U 3 k,- s (q^W2s(P -g- + Vl 2s {P - £,£>)] 

+fi>jks{q,i)[Ui>{P -q- £',£') + Ui(P - £,£')]} = M 2 f ijs (q), 



and 



mf + (gU + q 2 ±) 2 | Mj+(7i± | Vt + qh 



i - yi - yi 



yi 



2/2 



/jjfes (91, 92) 



(8) 



■£(-1/ 



V 1 + fijk 



p+ 



X { /^,- S (gi) [Vl 2s (P - £l - £2, £2) + ^2, (P - £1, £2)] 

(£l ) Pi {P - £1 - £2, £2) + C/i' (£ - £1 , £2)] 
+M,-s(£2)[^_* 2s (P - £1 - £2,£l) + V 2s (P - £2,£l)] 

+Uks (£2) [t/i (p - £1 - £2, £1) + Ui> (£ - £2, £1)] } + . . . 

= M 2 f ijks (q 1 ,q 2 ). 



We consider truncations of this system. 

A truncation to one boson leads to an analytically solvable problem. 7 
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The one-boson wave functions are 

P+ 



fij+(g) 



M 2 



mf+qi 



(9) 



fij-(g) = 



l-q+/P+ q+ /P+ 

Y,(-l) k+1 Zk)Ui(P-q,q) + J2(-l) k+1 ZkU k (P,q) 

k k 
P+ 

k 



M 2 



H9i 



■J2(-l) k+1 z k )V;(P-q,q). (10) 



l-q+/P+ q+/P+ 

Substitution into Eq. (6) yields 

(M 2 - m 2 )zi = g 2 nl{z - z{)J + g 2 m i (z () m - znn^Io 
+g 2 Ho[(zo ~ z\)mi + zairia - zim\)h, 



with 



In{M 2 



J(M 2 



dydq\ 



16tt 



E 



{mj/noY 



dydq 2 
16tt 



"j+gj _ ^1+9 i 2/(1 - y) n 
1-1/ y 



(_l)j+fc 



K 2 + <zi)/Mo 



M 2 



2/(1 -y) 2 



(ii) 



(12) 



(13) 



1-1/ y 

The presence of the PV regulators allows Iq and J to satisfy the identity 
^J(M 2 ) = M 2 I a (M 2 ). With M held fixed, the equations for z, can be 
viewed as an eigenvalue problem for g 2 . The solution is 

(M =p to )(M =p mi) 0i M =F m 



■9 = 



(14) 



(mi - m )(/U /i ± MJ )' z M^mi' 

An analysis of this solution is given in Ref. 7. 

In a truncation to two bosons, we obtain the following reduced equations 
for the one-boson-one-fermion wave functions: 8 

,2 



M 



1-2/ 



fijs(y,q±) 



(15) 



S 2 Iija(y,Q±) f i x 



2 /-l 

if^2 E / d y' d Q± J Siabs>(y><i->-'>y'> < i±)fabs>(y',q±) 

abs' " 



abs' 



+Y7^2E /. d y' d i'lJijlabs'(y^±;y' \q'±)fabs<(y' 
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where yfP+f lJ+ {q) = f ij+ {y,q±), \ f P+f lJ -{q) = fij-{y,q±)e i ' t ' , I is an 
analytically computable self-energy, and is a kernel determined by n- 
boson intermediate states. These reduced integral equations are converted 
to a matrix equation via quadrature in y' and q'±. The matrix is diagonal- 
ized to obtain g 2 as an eigenvalue and the discrete wave functions from the 
eigenvector. 

A useful set of quadrature schemes is based on Gauss-Legendre quadra- 
ture and particular variable transformations. The transformation for y' is 
motivated by the need for an accurate approximation to the integral J. 
This integral appears implicitly in the product of the Hamiltonian and the 
eigenfunction and is largely determined by contributions near the endpoints 
whenever the PV masses are large. The transformation for the transverse 
integral is chosen to reduce the range from infinite to finite, so that no 
momentum cutoff is needed. 

From the wave functions we can extract a structure function fssiy), 



f Bs (y)= J dq5(y-q+/P+) 
+ f f[dq IL J2S(y-q+/P+) 

J ^ _ 1 „ _ 1 



J2(-l) i+j f ijs (q) 



n—1 n—1 



^ ( _ 1)i+j+fe/ .. fes( ^ ) 

ijk 



(16) 



+ -.., (17) 



defined as the probability density for finding a boson with momentum frac- 
tion y while the constituent fermion has spin s. Typical results are plotted 
in Fig. 1. 



3. Feynman-gauge QED 

We apply these same techniques to QED. 9 The Feynman-gauge Lagrangian 
is 



i=0 



(18) 
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Figure 1. Bosonic structure functions in Yukawa theory, with a two-boson truncation 
{fs+'- solid; /s— : long dash) and a one-boson truncation (/s+: short dash; /s_: 
dotted). The constituent masses had the values mo = — 1.7//o, mi = Mi = 15/io- The 
resolutions used in the Gauss-Legendre method are K = 20 and AT = 30. 



where A» = ELo^f, ^ = ELo^, and if = d^A" - d v A£. The 
nondynamical fermion fields are constrained by 



i[(-iyd_^_-ieJ2A k _J2^ 

k j 

= -[(*7V)((-l) < ^i+-ieX;^fcxE^+)-(-l) < n»7Vi+]- 



(19) 



For the null combination = + V-'i-j this becomes 



id-tp- = — [(i7°7 )9j_V>+ — to7°?/> + ], 



(20) 



which is independent of A and can therefore be solved without inverting a 
covariant derivative. We then obtain the Hamiltonian without antifermion 
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terms as being 

P = E/ dp^^{-l)% s {p)hAp) (21) 

i,s 

i.fi 

+ E J d E d g{[ b UE)hAg)V^ 2s (p,q) 

+ bUp^hsigM^sibg)} Am- 2) + h - c ] . 

where e M = (—1, 1, 1, 1). The vertex functions U and V are given in Rcf. 9. 

The dressed electron state, without pair contributions and truncated to 
one photon, is 

|V) = 5>& i+ (P)|0> + E / dkfti s (k)bl(k)al(P- fc)|0), (22) 
with one-photon-one- electron wave functions 

f , m _ E.M^+^+fea 3 

a; 1— x 



1-x 



Substitution into P+P-|V>) = Af 2 |V>) yields 

m fe + ^1 — 2nik(mj + mi)x + rrijiriiX 



(M 2 - m 2 K = -1 / 2J(-l) i+fc+ % (25) 

(26) 



M 2 x(l - x) - m|(l - x) - ^ 2 :r - k\ 



This is the same form as in the one-boson Yukawa problem, with g 2 — > 
2e 2 and ii — » — 2Ji, and an analytic solution is again obtained. From 
this solution we can compute various quantities, including the anomalous 
magnetic moment. 9 
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4. A correlator in A/"=(2,2) SYM theory 

Reduction of Af=l SYM theory from four to two dimensions provides the 
action we need. In light-cone gauge (A_ = 0) it is 



= J dx + dx~ tr 



d+Xid-Xi + i6 T R d + e R + ie T L d-e L (27) 



+ l -{d-A+f +gA + J+ + V2gO T L e 2 (3 I [X I ,6 R ] + 9 ^\X u Xjf 

Here the trace is over color indices, the X[ are the scalar fields and the rem- 
nants of the transverse components of the four-dimensional gauge field , 
the two-component spinors Or and 8l are remnants of the right-moving 
and left-moving projections of the four-component spinor in the four- 
dimensional theory. We also define J+ = i[Xj,d-Xj] + 28^8 R , /3\ = a\, 
02 = 03, and e 2 = —102- 

The stress-energy correlation function for A/"=(8,8) SYM theory can be 
calculated on the string-theory side: 12 (T++(:r)T++(0)) = (N? /2 /g) X - 5 . 
We find numerically that this is almost true in W=(2,2) SYM theory 14 

To compute the correlator, 13 we fix the total momentum P+, compute 
the Fourier transform, and express the transform in a spectral decomposed 
form 

F(P+,x~) - ±.(T++(P+,x+)T++(-P+, 0)) (28) 

= E^(°i T++ ( p+ )i i ) e " Pra;+ ( i i T++ (- p+ ' )i )- 

i 

The position-space form is recovered by Fourier transforming with respect 
to the discrete momentum P + = Ktt/L, where K is the integer resolution 
and L the length scale of DLCQ. 2 This yields 

I T 2 / T + \ 2 /If 4 

F(x-,x + ) = J2\-W ++ (K)\i) [—) ^^K i (M l V2-x+x^). (29) 



We then continue to Euclidean space by taking r = V2x+aP to be real. The 
matrix element (L/ir)(0\T ++ (K)\i) is independent of L. Its form can be 
substituted directly to give an explicit expression for the two-point function. 
The correlator behaves like 1/r 4 at small r: 



F( X -, X +)= N °^ 2 + n '\ l-l,K). (30) 
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For arbitrary r, it can be obtained numerically by either computing the 
entire spectrum (for "small" matrices) or using Lanczos iterations (for 
large). 13 

In Fig. 2, we plot the log derivative of the scaled correlation function 14 



At small r, the results for / match the expected (1 — 1/K) behavior. 



Figure 2. The log derivative of the scaled correlation function / defined in Eq. (31) of 
the text. The resolution K ranges from 3 to 12. For even K, f becomes constant at 
large r and the derivative goes to zero. 

At large r the behavior is different between odd and even K, but as K 
increases, the differing behavior is pushed to larger r. For even K, there 
is exactly one massless state that contributes to the correlator, while there 
is no massless state for odd K. The lowest massive state dominates for 
odd K at large r; however, this state becomes massless as K — > oo. In 
the intermediate-r region, the correlator behaves like r -4 ' 75 , or almost r~ 5 . 
The size of this intermediate region increases as K gets larger. 

5. AT— (1,1) SYM theory at finite temperature 

In this case, the Lagrangian is L = Tr (-\F llv F> lv + i^j^D^), with 
Fy, v = d^Av - dvAy, + ig[A^,A v ] and = + ig[A^]. The super- 
charge in light-cone gauge is Q~ = 2 3//4 <? J dx~ (i[<j>, d-(f>] + 2i[>ip) dZ 1 i> ■ 
From the discrete form we can compute the spectrum, which at large- N c 




(31) 
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represents a collection of noninteracting modes. With a sum over these 
modes, we can construct the free energy at finite temperature from the 
partition function 17,15 e~ p °/ T . 

The one-dimensional bosonic free energy is 

and the fermionic free energy is 

= -— E f ^°-rP=m ln ( x + e " 0/T ) • (33) 



The contributions from the K — 1 massless states in each sector are 

*S = -^^^VT 9 , ^ = - { ^^VT\ (34) 

The total free energy, with the logs expanded as sums and the po integral 
already performed, is 

rW.JJ^rt-^Xf^KMtm. ,35, 

n=l 1=0 

The sum over I is well approximated by the first few terms. We can 
represent the sum over n as an integral over a density of states: ^2 n — > 
/ p{M)dM and approximate p by a continuous function. The integral over 
M can then be computed by standard numerical techniques. We obtain p by 
a fit to the computed spectrum of the theory and find p{M) <~ exp(M/Tn), 
with Th ~ 0.845 \A7g 2 /Vc, the Hagcdorn temperature. 18 From the free 
energy we can compute various other thermodynamic functions up to this 
temperature. 15 



6. Future work 

Given the success obtained to date, these techniques are well worth con- 
tinued exploration. In Yukawa theory, we plan to consider the two-fermion 
sector, in order to study true bound states. For QED the next step will be 
inclusion of two-photon states in the calculation of the anomalous moment. 
For SYM theories, we are now able to reach much higher resolutions, by 
computing on clusters. This will permit continued reexamination of theories 
where previous calculations were hampered by low resolution, particularly 
in more dimensions. Earlier work on inclusion of fundamental matter 19 can 
be extended to three dimensions and modified to include finite- N c effects, 
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such as baryons with a finite number of partons and the mixing of mesons 
and glueballs. For all of this work, the ultimate goal is, of course, the 
development of techniques sufficient to solve quantum chromodynmics. 
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